Study of glass transition kinetics of As2S3 and As2Se3 by ultrafast differential scanning calorimetry
Zhang Fan1, Chen Yimin1, 2, †, Wang Rongping1, Shen Xiang1, Wang Junqiang3, 4, , Xu Tiefeng1
Laboratory of Infrared Material and Devices & Key Laboratory of Photoelectric Materials and Devices of Zhejiang Province, Advanced Technology Research Institute, Ningbo University, Ningbo 315211, China
Department of Microelectronic Science and Engineering, Faculty of Science, Ningbo University, Ningbo 315211, China
CAS Key Laboratory of Magnetic Materials and Devices & Zhejiang Province Key Laboratory of Magnetic Materials and Application Technology, Ningbo Institute of Materials Technology & Engineering, Chinese Academy of Sciences (CAS), Ningbo 315201, China
Center of Materials Science and Optoelectronics Engineering, University of Chinese Academy of Sciences, Beijing 100049, China

 

† Corresponding author. E-mail: chenyimin@nbu.edu.cn jqwang@nimte.ac.cn

Abstract

Ultrafast differential scanning calorimetry (DSC) was employed to investigate the glass transition kinetics of As2S3 and As2Se3. By using the Arrhenius method, a fragility index of ∼22 can be estimated in both As2S3 and As2Se3. However, when the scanning rate is more than , non-Arrhenius behavior can be observed in such “strong” liquids where the Vogel–Fulcher method is more accurate to describe the glass transition kinetics. The fragilities of As2S3 and As2Se3 glasses are thus extrapolated as 28.3 ± 1.94 and 23.7 ± 1.80, respectively. This indicates that, As2Se3 glass has a better structural stability and it is a better candidate for device applications.

1. Introduction

Structural relaxation is intrinsic in amorphous materials which results in time-dependent physical properties or physical ageing.[1] When a glass is quenched from liquid state, it is characterized by a lack of thermodynamic equilibrium. Therefore, once the glass is maintained at a temperature below its glass transition temperature (Tg) in a period of time, its properties are changed more or less, implying the physical ageing is closely correlated to structural stability. An effective way to evaluate the structural stability of glass forming liquids is the conception of “fragility” introduced by Angell.[2] The fragility index, which is determined by the change of relaxation time or viscosity at the glass transition temperature, has been proven useful in the study of structural relaxation and glass transition kinetics in the glass forming liquids.[3]

Some measurement methods have been proposed to obtain the quantitative fragility index, such as relaxation time[4] and viscosity.[5,6] However, these experimental processes are difficult to control and time-consuming. Another effective way to estimate the fragility is equation parameterizations via glass transition kinetics study. Conventional differential scanning calorimetric (CDSC) investigation on the glass transition kinetics suggests that, the scanning rate-dependent Tg measurement during heating or cooling obeys an exponential law.[3] Nevertheless, this is an imprecise method due to the narrow range of the scanning rate. Recently, a novel ultrafast differential scanning calorimetry (UDSC) has been developed and successfully used for the study of the glass transition kinetics in Te-based chalcogenides, such as Ge2Sb2Te5,[7] Ag–In–Sb2Te,[8] and Ge–Te.[9,10] The scanning rate can be extended up to during the measurements, thus, it is possible to investigate the abnormal glass transition kinetics hidden in some supercooled liquids and estimate the liquid fragility more accurately using UDSC.

In present paper, we employed UDSC to investigate the glass transition kinetics of As2S3 and As2Se3. We precisely measured the fragility indexes of S- and Se-based arsenic glasses using the Vogel–Fulcher method. It was found that, As2Se3 has smaller fragility index with a better structural stability.

2. Experimental methods
2.1. Sample fabrication

As2S3 and As2Se3 glasses were fabricated by the melt-quenching method. The high purity raw materials of As (5 N), S (6 N), and Se (5 N) were weighted, and about 20-g mixtures were loaded into quartz ampoules (ϕ = 10 mm), which were subsequently sealed under a vacuum less than 10−3 Pa. The ampoules were shifted into a rocking furnace and heated slowly to around 700 °C for 12 h, then homogenized at 500 °C for 0.5 h. The ampoules were quenched in water to allow glass formation. Afterward, the vitreous samples were annealed at 30 °C below their glass transition temperatures for 12 h to relax the internal mechanical stress. Then the glass rods were cut and polished for further measurements.

2.2. Calorimetric measurements

CDSC measurements were carried out in TA Q2000 with a temperature accuracy of ±1 K. In each CDSC measurement, about 10-mg sample was enclosed into aluminum pan and heated at a scanning rate of under N2 atmosphere. UDSC measurements were performed in Mettler Toledo Flash DSC 1. In each UDSC measurement, a tiny sample with a mass of hundreds of nanogram was put on a chip sensor, and then heated at different heating rate from to under an Ar atmosphere.

2.3. Thermal lag in UDSC measurement

In this work, we concentrate on the glass transition kinetics investigated by UDSC. Thermal lag is a major concern since the heating and cooling rates are ultrafast in UDSC. Biot number is usually used to estimate the potential thermal lag at the heating process in the studies. The Biot number is defined as, ,[7] where L is the sample thickness, h is the heat-transfer coefficient between the heater and surface (a typical value of is employed in here),[11] and κ is the thermal conductivity of the sample. The value of κ at room temperature is in a range from to for As2S3,[12] and more than for As2Se3.[13] It should be noted that, we used the crushed flakes with a thickness of 2–5 microns for UDSC measurements. Thus, the Biot number can be estimated as 0.028–0.07. Previous investigations indicating that the thermal lag is negligible in UDSC measurement if the Biot number is less than 0.1.[7,9] Therefore in the present paper, the effect of thermal lag is negligible.

3. Results and discussion

Figure 1(a) and 1(b) show the subtracted UDSC traces for As2S3 glass with different heating rates from to and to , respectively. A peak-like endothermic stage, which consists of an endothermic stage and the excess heat capacity peak determined by the atom translational motions in glass transition,[14] can be found in each subtracted trace. The onset glass transition temperature (Tg) increases with increasing heating rate. A similar scenario was found in As2Se3 glasses. Figure 2 depicts the conventional DSC traces of As2S3 and As2Se3 glasses obtained at a heating rate of . Tg is 481 K and 467 K for As2S3 and As2Se3, respectively. Those values are in excellent agreement with the previous data reported.[1518]

Fig. 1. Subtracted UDSC traces for As2S3 with the heating rate are (a) to and (b) to . The subtracted UDSC traces for As2Se3 glasses with the heating rate range from (c) to and (d) to .
Fig. 2. Conventional DSC traces tested at a heating rate of for As2S3 and As2Se3.

T g obtained at a large range of the heating rate makes it possible to investigate the glass transition kinetics. Figure 3 describes the relationship between natural logarithm of the heating rate and reciprocal Tg (1000/Tg). It shows an Arrhenius behavior, especially at the heating rate less than in both As2S3 and As2Se3. However, when the heating rate is higher than (the data above divergent point in Fig. 3), a non-Arrhenius behavior is obviously observed in As2S3, but slightly presented in As2Se3. We first used a linear fit method to estimate the glass transition kinetics at a low heating rate range (light yellow shade in Fig. 3). As shown in Fig. 3, the dashed lines with the slopes ( , where E a is the activation energy and R is the gas constant) of −24.916×103 and −23.505×103 can be obtained for As2S3 and As2Se3, respectively. According to the relationship, ,[19] where Tg is the glass transition temperature, i.e., 481 K and 467 K for As2S3 and As2Se3, their fragilities can be estimated as 22.5 and 21.9. Thus, As2S3 and As2Se3 glasses possess similar fragility index that is close to 22, and they are the strong liquids that can remain thermodynamically stable over a wide temperature range.[20] 4

Fig. 3. The natural logarithm of heating rate, ln ϕ, is plotted against 1000/Tg for (a) As2S3 and (b) As2Se3. The black dashed lines are the results of linear fitting, and the slope of the fitted line is . Green triangle and purple diamond in Fig. 3(a) represent the As2S3 data from Espeau et al.,[15] and Street and Yoffe.[16] Orange triangle and olive diamond in Fig. 3(b) represent the As2Se3 data from Blachnik and Hoppe,[17] and Chen et al.[18]
Fig. 4. The relation between Tg and heating rate of (a) As2S3 and (b) As2Se3. The dashed curves were fitted by Vogel–Fulcher equation.

Moving to the deep yellow shade in Fig. 3, the linear fit cannot well depict the relation between ln ϕ and 1000/T g when ϕ is more than . Such non-Arrhenius behavior observed in the ultrafast DSC measurements enables the Vogel–Fulcher method to be used for calculating the fragility index more accurately. The relation between Tg and ln ϕ in the Vogel–Fulcher method is,[21]

where B is the time scale parameter in glass forming system, D is the strength parameter represents how close the liquid obeys Arrhenius law, and T0 is the asymptotic value of Tg within the limit of infinitely slow heating (or cooling) rate.[22] Afterward, the fragility index can be estimated as,[4]
where T g is glass transition temperature obtained at a heating rate of (see Fig. 2). The obtained parameters were listed in table 1. Apparently, As2Se3 has a large D index of 54.9. This indicates As2Se3 is the stronger liquid and obeys the Arrhenius law. Taking the parameters D, T0, and Tg into Eq. (2), the fragility m can be calculated as 28.3 ± 1.92 and 23.7 ± 1.80 for As2S3 and As2Se3, respectively. Compared with the results estimated by linear fit, both fragilities become larger.

Table 1.

The parameters fitted by Vogel–Fulcher equation and fragilities for As2S3 and As2Se3.

.

Using the data of viscosity close to Tg, the fragility index of the As2Se3 liquid was reported by Málek et al. as .[23] Brazhkin et al. reported that, the viscosities are 105, 10, 0.1, and 0.01 Pa s at 350, 500, 700, and 1000 °C, respectively, for As2S3 liquid.[24] Combining with the Mauro–Yue–Ellison–Gupta–Allan (MYEGA) viscosity model,[6] a rough value of fragility was estimated as ∼46 for As2S3. It can be seen that, the fragility index of As2Se3 is smaller than that of As2S3, which is in line with our result in this work. In addition, the fragility indexes estimated from viscosity change for both As2S3 and As2Se3 are all larger than those estimated from glass transition kinetics. Depending on the viscosity data, these arsenic glasses may consider as the materials with intermediate kinetics, which have the fragility index between the strong and fragile pattern, rather than the “strong” materials. However, we know the strong glass forming liquids are characterized by covalent directional bonds that form a spatial network, such as the SiO2.[23] Recently, by using the conventional DSC, the fragility of As2S3 was roughly estimated as 27,[25] which is very close to our result. Thus, the typical covalent glasses, i.e., As2S3 and As2Se3, can be considered as the “strong” glass forming liquids.

4. Conclusion

We have studied the glass transition kinetics of As2S3 and As2Se3 glasses using novel UDSC. According to the linear fitting results, it yields the fragility index of 22.5 and 21.9 for As2S3 and As2Se3, respectively. However, the non-Arrhenius behavior was observed when the heating rate is more than , indicating that the linear fit is inappropriate in the large heating rate. Thus, the Vogel–Fulcher method was employed to estimate the fragility index more precisely. It was found that, the fragilities for As2S3 and As2Se3 are 28.3 ± 1.94 and 23.7 ± 1.80, respectively, which are all larger than those estimated from the linear fit. This fact illustrates that, even for the “strong” glasses, large range of scanning rate (such as those used in UDSC) is helpful to estimate fragility index more accurately for glass transition kinetics.

Reference
[1] Qiao J C Wang Q Crespo D Yang Y Pelletier J M 2017 Chin. Phys. 26 016402
[2] Angell C A 1995 Science 267 1924
[3] Wang L M Velikov V Angell C A 2002 J. Chem. Phys. 117 10184
[4] Böhmer R Ngai K L Angell C A Plazek D J 1993 J. Chem. Phys. 99 4201
[5] Evenson Z Schmitt T Nicola M Gallino I Busch R 2012 Acta Mater. 60 4712
[6] Mauro J C Yue Y Ellison A J Gupta P K Allan D C 2009 Proc. Natl. Acad. Sci. USA 106 19780
[7] Orava J Greer A L Gholipour B Hewak D W Smith C E 2012 Nat. Mater. 11 279
[8] Orava J Hewak D W Greer A L 2015 Adv. Function. Mater. 25 4851
[9] Chen Y Wang G Song L Shen X Wang J Huo J Wang R Xu T Dai S Nie Q 2017 Crystal Growth & Design 17 3687
[10] Chen Y Pan H Mu S Wang G Wang R Shen X Wang J Dai S Xu T 2019 Acta Mater. 164 473
[11] Şahin H M Kocatepe K Kayıkcı R Akar N 2006 Energy Convers. & Manag. 47 19
[12] Kim S W Yu H Kang C H Lee S H Kim J C 2004 Int. J. Thermophys. 25 1505
[13] González-Romero R L Antonelli A Chaves A S Meléndez J J 2018 Phys. Chem. Chem. Phys. 20 1809
[14] Ke H B Wen P Wang W H 2012 AIP Adv. 2 041404
[15] Espeau P Tamarit J L Barrio M López D Perrin M Allouchi H Céolin R 2006 Chem. Mater. 18 3821
[16] Street R A Yoffe A D 1972 J. Non-Cryst. Solids 8�?0 745
[17] Blachnik R Hoppe A 1979 J. Non-Cryst Solids 34 191
[18] Chen P Holbrook C Boolchand P Georgiev D Jackson K Micoulaut M 2008 Phys. Rev. 78 224208
[19] Chen Y Mu S Wang G Shen X Wang J Dai S Xu T Nie Q Wang R 2017 Appl. Phys. Express 10 105601
[20] Bi Q L Y J 2014 Chin. Phys. Lett. 31 106401
[21] Zhang B Wang R Zhao D Pan M Wang W 2004 Phys. Rev. 70 224208
[22] Brüning R Samwer K 1992 Phys. Rev. 46 11318
[23] Málek J Shánělová J 2005 J. Non-Cryst. Solids 351 3458
[24] Brazhkin V V Kanzaki M Funakoshi K Katayama Y 2009 Phys. Rev. Lett. 102 115901
[25] Qiao A Wang P Tao H Yue Y 2018 J. Non-Cryst. Solids 500 225